Root system architecture of historical spring wheat cultivars is associated with alleles and transcripts of major functional genes

We evaluated root system architecture (RSA) of a set of 58 historical spring wheat cultivars from Pakistan representing 105 years of selection breeding. The evaluations were carried out under control and water-limited conditions using a high-throughput phenotyping system coupled with RhizoVision Explorer software. The cultivars were classified into three groups based on release year as cultivars released pre-1965, released between 1965 and 2000, and cultivars released post-2000. Under water-limited conditions a decline in 20 out of 25 RSA component traits was observed in pre-1965 cultivars group. Whereas cultivars released after the 1965, so-called green revolution period, showed a decline in 17 traits with significant increments in root length, depth, and steep angle frequency which are important root traits for resource-uptake under water-limited conditions. Similarly, cultivars released after 2000 indicated an increase in the number of roots, depth, diameter, surface area, and steep angle frequency. The coefficient of correlation analysis showed a positive correlation between root depth and yield-related traits under water-limited conditions. We also investigated the effects of green-revolution genes (Rht1) and some phenology-related genes such as DRO1, TaMOR, TaLTPs, TaSus-2B on RSA and identified significant associations of these genes with important root traits. There was strong selection pressure on DRO1 gene in cultivated wheat indicating the allele fixed in modern wheat cultivars is different from landraces. The expression of DRO1, and TaMOR were retrieved from an RNAseq experiment, and results were validated using qRT-PCR. The highest expression of DRO1 and TaMOR was found in Chakwal-50, a rainfed cultivar released in 2008, and MaxiPak-65 released in 1965. We conclude that there is a positive historic change in RSA after 1965 that might be attributed to genetic factors associated with favored RSA traits. Furthermore, we suggest root depth and steep angle as promising traits to withstand water-limited environments and may have implications in selection for breeding. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-03937-7.


Introduction
One of humanity's greatest challenges in 21 century is feeding the ~ 9 billion people by 2050 with continuously decreasing natural resources [1]. Extreme weather events associated with climate change have threatened crop production patterns across the globe [2]. Drought and flooding events in response to climate extremes may increase, and it is critical to devise possible strategies to enhance crop production in the forthcoming decades [3,4]. Wheat (Triticum aestivum L.) is one of the major staple food crops worldwide supplying 20% calories to the global population. Wheat is sensitive to drought and heat stress, that are critical environmental adversity affecting final crop yields [5]. Breeding for drought tolerant genotypes with agronomic and adaptive traits is crucial for increasing productivity and food security among wheatgrowing communities.
Root system architecture (RSA) is important for water and nutrient uptake facilitating photosynthesis and improving grain yield. It represents a fair benchmark for manipulation under soils poor in nutrition improving crop productivity [6,7]. RSA refers to the spatial organization of root structure that includes the root's number, length, tip number, emergence angles, depth, width, convex hull area, and root mass center [8]. To optimize the nutrient and water uptake, understanding and selection of unique RSA traits and identification of their underlying genes are equally important as above-ground components. There are specific root architectural traits that benefit crop yield by improving soil resource acquisition. For instance, a deeper RSA can extract more water and nutrients under drought conditions relative to a wider RSA in optimal conditions. Root traits contributing to plant productivity under drought stress such as long specific root length, small fine root diameter, and high root length density have also gained significant attention [9]. In wheat, seminal roots (lateral roots originating from radicle) largely determine the root architecture at the adult plant stage. Several traits of these roots such as root growth angle, seminal root number, and length can be conveniently assessed at early growth stages [10].
Recently, a genetic framework of 38 genes underpinning RSA in cereals was reported [11]. For instance, the reduced height (Rht) genes responsible for the green-revolution had a significant impact on root traits and coleoptile length [12,13]. Overexpression of TaMOR (More roots in wheat) in rice plants resulted in longer main panicle, more crown roots, a higher number of primary branches and an increased grain yield. Root depth is an important trait allowing better access to nutrients and water stored in deeper soil layers thus enhancing yield [14]. DRO1 influences root growth angle by modulating root gravitropic response. The deep rooting allele has been found to increase grain yield under drought stress [15]. Three other major QTL namely DRO2, DRO3, and qSOR1 have been reported to control RSA under waterdeficient conditions [14]. Rambla et al. [16] reported a single plant selection (SPS) approach for introgression of root traits into elite wheat germplasm. The approach combines phenotypic selection (root angle and biomass), marker assisted selection (MAS) using KASP markers for qRDM-5B, a root biomass QTL, and speed breeding to accelerate the breeding cycles.
The main challenge in breeding programs is the lack of high-throughput phenotyping platforms for root traits that can offer a proxy for field performance. Despite that, genetic improvement in RSA is an underappreciated route to a more efficient and productive wheat crop. Evidence suggests that historical improvements in wheat productivity are linked to changes in RSA [17]. Hence, designing a root system tailored to target the edaphic environment via modification of genes underpinning favorable root traits represents an ideal breeding strategy for direct selection [4]. The development and release of new adaptable varieties with desired root systems could be a promising strategy to surmount unfavorable environmental conditions. Given the above factors, this study was conducted to investigate phenotypic variability in RSA traits of historical wheat cultivars representing the breeding progress since 1911 to 2016 in Pakistan and to identify the associations of unique RSA with phenology, yield-related traits, and allelic variations and expression of important functional genes.

Germplasm
A panel of 58 bread wheat cultivars released in Pakistan (1911Pakistan ( -2016 was selected for this study. The cultivar name, year of release, and pedigree are given in Table 1. A subset of this collection was selected for RNA-sequencing and qRT-PCR validation. The schematic overview of the experiments is given in Fig. 1. All the experimentation was conducted on wheat plants, so no formal identification of plants was required, and no voucher specimen was needed to deposit in herbarium. The seeds of the cultivars used in this experiment can be accessed from National Genebank, Plant Genetic Resources Institute, National Agriculture Research Center (NARC), Islamabad, Pakistan.

Phenotyping experiments Evaluation of drought stress response in cigar rolls
Seeds were surface sterilized with 2% NaOCl (BioChem; active chlorine 4.0-6.0) and subsequently rinsed with distilled water three times. Ten seeds of each cultivar were placed on a Whatman quantitative ashless filter paper (Grade No.42) fixed in Petri plate and kept in dark for germination. After the emergence of coleoptiles 5 days after sowing (DAS), five uniform seedlings were transferred to germination papers rolled in cigars configuration. The cigar rolls were then placed in 1 L glass beakers totaling eight cigars in each beaker filled with 300 mL distilled water and kept as control. For drought stress treatment, 150 mL distilled water was maintained every day in   [18]. Seedlings were removed 2 weeks after treatment, roots were separated and imaged using Rhizo-Vision Crown hardware platform coupled with a machine vision camera. RhizoVision crown hardware [19] controlled by Rhizo-Vision Imager software was used for root image acquisition. The hardware consisted of a backlit imaging box with 65 cm × 66 cm × 91 cm dimensions that produced near binary images. The imaging box was constructed from T-slotted aluminum profiles. LED edge-lit flat panel light (65 cm × 65 cm) affixed at the back of the box provided a white background while roots gave the foreground near-black color. A root holder (22 cm × 30 cm) was constructed at the top of the box and a door handle was attached on the top of the root holder panel to assist in the placement and removal of roots. A monochrome machine vision camera (Basler acA5472-17 μm USB 3.0) with a CMOS sensor (Sony IMX183) was used to capture images. For root image analysis, an open-source software RhizoVision Explorer was used [20].

Field experiments
The same set of 58 wheat cultivars was evaluated for two-year (2018 and 2019) in the field at National Agriculture Research Centre (NARC), Islamabad, Pakistan using a randomized complete block design (RCBD) with two replications. The NARC site is located at 33°43′N 73°04′E and has soil electrical conductivity of 0.24dS/m. The plot size of 6 rows of 6 m length with 16 cm rowrow distance was maintained. Phenotyping data for various morphological parameters including tillers per plant (TPP), plant height (PH) in cm, spike length (SL) in cm, spikelet per spike (SNPS), grains per spike (GPS), thousand grain weight (TGW), grain length (GL), grain width (GW), grain diameter (GD), and grain yield (GY) in kg per m 2 per plot were recorded as described in our previous study [21]. The experiments were conducted under standard institutional policies and no biosafety permission was needed for this experimentation. No permission was needed to collect plants and plant materials. All procedures were conducted in accordance with the institutional guidelines.

Genotyping for functional genes
Genomic DNA of cultivars was extracted from leaf tissues using the phenol-chloroform method. The KASP markers for genes Rht-B1, Rht-D1, and TaSus2-2B, were used from our previous work [22]. The KASP assay was prepared using 2 μl of 50 ng/μl DNA as template, 2.5 μl of 2X KASP master mix, 0.07 μl of KASP assay mix, and 2.5 μl of PCR H 2 O. The PCR was performed in 384-well formats (S1000, Thermal Cycler, USA) following conditions described by [21].

RNAseq and qPCR-based expression of DRO1 and TaMOR
A subset of 12 wheat cultivars was subjected to RNAseq to identify the genotypic variation in gene expression in seedling leaf and roots (PRJNA863398). For this purpose, seedling leaf and root tissues in triplicate were collected and subjected to total RNA extraction. Briefly, RNA was extracted using EasyPure Plant RNA Kit (ER301-01) following manufacturer instructions. RNA quality was checked on 1% Agarose gel and quantified using Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). The samples were sent to Beijing Genomics Institute (BGI), China for sequencing. In summary,mRNAs were isolated and fragmented from total RNA using the oligo (dT) method for cDNA synthesis. The 50-bp single-end sequencing libraries were constructed, and sequencing was performed on BGISEQ-500 platformusing standard protocols. Adapters with unknown bases (N's > 5%) and low quality were removed from raw reads to produce 'clean data' as FastQ data files using a quality control software, SOAPnuke version 2.1.6.

Identification of differentially expressed genes
High quality single-end reads were mapped to the bread wheat reference genome (IWGSC, INSDC Assembly GCA_900519105.1) using HISAT2 (Hierarchical Indexing for Spliced Alignment of Transcripts) software version 2.2.1 [23]. Alignment of the reference sequence with reads were performed using Bowtie [24]. Quantification of the reads was performed using featureCounts software program. Differentially expressed genes (DEGs) were identified using DeSEQ2 in R version 4.1.1. DEGs were then filtered based on their adjusted p-value. The threshold value for filtering was set at 0.1. To analyze the expression of DRO1 and TaMOR genes in wheat genotypes, the heatmaps were constructed using pheatmap package in R version 4.1.1.

Validation of DEGs using qRT-PCR
For expression analysis of DRO1 and TaMOR genes, a subset of 12 cultivars selected based on high and least root length was germinated in cigar rolls. A total of five surface-sterilized seeds of each cultivar were placed between two germination papers and rolled in a cigar configuration. Two batches of fifteen cigar rolls totaling ten plants per genotype were placed in beakers containing Hoagland's nutrient solution. The experiment comprised of two treatments i.e., well-watered (WW) and water-limited (WL) conditions. For well-watered treatment, 200 mL of Hoagland's solution was added while 100 mL nutrient solution was supplied to maintain WL conditions. The nutrient solution was supplied daily to maintain a volume of 200 mL and 100 mL for WW and WL conditions respectively. The roots were harvested 15 days after germination for RNA extraction from each treatment.

RNA isolation, cDNA synthesis, and qRT-PCR analysis
Total RNA was extracted from a set of 12 wheat cultivars using EasyPure Plant RNA Kit (ER301-01) following manufacturer instructions. RNA quality was checked on 1% Agarose gel and quantified using Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). RNA samples with 1.8-2.1 A260/A280 values were selected for cDNA synthesis. Genomic DNA was removed, and cDNA was synthesized using ABClonal ABScript III RT Master Mix with gDNA remover. The reverse transcription reaction system consisted of 4 μl 5X ABScript III RT Mix, 1 μl 20X gDNA remover mix, 1 pg-1 μg of total RNA, and 13 μl nuclease-free water making up to 20 μl. The conditions for the reaction were 37 °C for 2 minutes, 55 °C for 15 minutes, 85 °C for 5 minutes, and 4 °C to hold. The products were quantified using Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA) and stored at − 20 °C for subsequent qRT-PCR reaction.
The coding sequences of DRO1 and TaMOR were retrieved from NCBI and were used to design common primers for all three homoeologues using Primer-BLAST (https:// www. ncbi. nlm. nih. gov/ tools/ primer-blast).
A previously designed primer pair for TaActin was used as an internal control. The sequences of primers used are given in Table 2. The transcription levels of DRO1 and TaMOR were quantified using Livak method in a CFX384 Real-Time detection system (Bio-Rad). The reaction components were as follows: 10 μl 2X Universal SYBR Green Fast qPCR, 1 μl cDNA product (40 ng/ μl), 0.4 μl forward primer (10 μM), 0.4 μl reverse primer (10 μM) and 8.2 μl nuclease-free water to make a final volume of 20 μl. The two-step reaction conditions were as follows: 1 cycle at 95 °C for 3 minutes and cycles at 95 °C for 5 seconds, 60 °C for 30 seconds. The amplification and melting curves were confirmed after the reaction and then a standard curve for quantitative analysis was made.

Statistical analysis
Data were analyzed using a generalized linear mixed model. Mean comparisons were made with Tukey's test at a p = 0.05 significance level. All statistical analyses were performed using Jamovi version 1.8. PCA analysis of root and yield related traits was performed using two principal components in Jamovi version 1.8. Correlations among root and yield traits were determined using Pearson correlation coefficient using GGally package in R 4.1.1. The student's t-test was used to investigate the significant association of a marker with the target trait at a threshold probability p < 0.05 [25].

Trends and genotypic differences in root traits over time
Wheat RSA traits were phenotyped in the Rhizovision Hardware Crown platform using the Cigar roll method. The set of fifty-eight wheat cultivars was divided into three groups based on breeding eras as pre-1965 (n = 6), 1965-2000 (n = 19), and post-2000 (n = 33). Descriptive statistics of each group with mean and range of all root traits under controlled and water-limited condition is given in Table 3 Significant variations (P < 0.05) among genotypes were found for all root traits except maximum diameter (MaxD) and width to depth ratio (WDR). The mean squares of each trait are given in Table 4. For treatment and variety x treatment interaction, significant differences were also found in 12 and 14 root traits, respectively. The PCA biplots exhibiting relationships among genotypes and various root and yield-related traits are illustrated in Fig. 3. The first two axes together explained 42.5% variation in control and 42.8% variation under drought stress. The angle of the vectors indicated marked differences among the association of various root traits with yield parameters in control and drought conditions. The SAF and WDR indicated remarked differences under control and water-limited conditions (Fig. 3A, B).

Allelic variations for functional genes and their association with root traits
The KASP markers were used to identify the allelic variations for 19 functional genes. Allelic effects and frequencies of some important genes are given in Table 5. Among Rht1 genes, Rht-D1a had the highest frequency (89.6%). Among TaSus2 The associations of the alleles with root traits demonstrated remarked differences (Fig. 4). Rht-B1 and TaMOC had no significant effect on any root trait tested. The TaSus2-2B, TaSnRK2.3-B1, TaSnRK2.9-5A, TaDreb-B1, 1fehw3, and TaPPH were significantly associated with root depth under water-limited conditions. TaDreb-B1 and 1fehw3 also had significant association with root volume. The TaSnRK2.3-B1 and TaLTPs were significantly associated MaxW. None of the alleles were found associated with root traits related to number and diameter. Individual Rht-B1 and Rht-D1 did not show any association with any traits except when combined in haplotypes. Rht1 haplotypes had significant effects on TRL, MaxW, NtA, SA, VD, convex area (CA), average diameter (AD), perimeter (P) projected area diameter (PAD), and surface area diameter (SAD) under control. No significant association of any allele was found with D under control condition.

Analysis of genetic variations at major root biomass QTL qRDM-5B
A major QTL, qRDM-5B, for dry root biomass was identified on chr5B [26]. Single plant selection protocol including the haplotype of this QTL was suggested to be used for selection of desired RSA in wheat [16]. Since the 660 K data was available for this collection, all the SNPs within the 653.8 to 654.9 Mb were retrieved and haplotypes were constructed based on linkage disequilibrium. Table 4 Analysis of Variance (ANOVA) for root traits of wheat cultivars analyzed under water-limited conditions * Significant (P < 0.05); ** Significant (P < 0.01); *** Significant (P < 0.001); ns Non-significant (P > 0.05); df: degree of freedom; All the values are mean squares In total, 21 different haplotypes were identified. Out of these, chr5B-hap20 consisted of 9 SNPs and four different haplotype alleles (Fig. 5a). The frequency of two alleles was very rare and were present only in two cultivars and were excluded from the association analysis. Two haplotype alleles of chr5B-hap-20 viz. CAG TTG TA  Table 2 (n = 13) and TGA CCA CGT (n = 26) were used for the association analysis. The student's t-test was significant where the former haplotype allele was associated with higher root dry weight (Fig. 5b), and the latter haplotype allele was associated with lower root dry weight (Fig. 5c).

Nucleotide variations and gene expression analysis of DRO1 and TaMOR
DRO1 and TaMOR genes were selected for expression analysis in roots. A 660 K SNP array data on this cultivars collection was used to extract SNPs within DRO1 and TaMOR genes ( Table 6). In DRO1, a total of 6 SNPs were identified out of which AX-109484887, and AX-108936574 caused missense mutation on DRO1-5B. Since no missense mutation was identified in DRO-5A, we then used SnpHub portal to identify any missense mutation in global wheat collection based on exome capture data. A SNP was identified with a very rare missense mutation in global wheat collection, and geographic map of DRO-A1 haplotypes is shown as Fig. 6a. Knetminer   Table 2 gene network of DRO1 showing associated traits and SNPs in Cadenza TILLIG population is shown (Fig. 6b). Since all the SNPs were fixed in the 58 cultivars panel with no allelic variation, therefore 660 K SNP array data was retrieved for 159 cultivars and landraces from Pakistan for in-depth analysis of SNP states. It was revealed that all the SNPs in DRO1 were fixed in wheat cultivars (113 out of 159) and alternate allele was identified in landraces (46 out of 159), and some older cultivars like Chenab-70, Local White and Kharchia. While AX-110640724 was completely fixed in all 159 accessions with no variation. Four SNPs were identified in TaMOR, of which AX-108991022 was present on both A-and Dhomeologue caused missense variations, and only three landraces had the nucleotide change.
Transcriptome analysis of twelve wheat genotypes showed significant differences in expression patterns of genes (Fig. 7a). DRO1 was expressed in both leaf and roots with variations in leaf tissues among genotypes. In roots, the highest expression of DRO1 in all genotypes was found in B homoeologue (TraesCS5B02G210500) with negative expression in A (TraesCS5A02G213300)  and D homoeologues (TraesCS5D02G218700) relative to B homoeologue. TaMOR exhibited no expression in leaf tissues in all genotypes, hence, differentially expressed in root tissues. The highest expression in A homoeologue (TraesCS4A02G415400) was found in Chakwal-50 cultivar with 0.70 tpm value whereas in B homoeologue (TraesCS4B02G316200) was found in C-518 (0.88 tpm) following Zincol-16 (0.45 tpm). The highest expression in D homoeologues (TraesC-S4D02G312800) was found in Pak-81 (1.10 tpm) following GA-2002 (0.54 tpm). qRT-PCR for DRO1 expression (Fig. 7b) indicated the highest normalized expression ratio in Chakwal-50 (29.45) whereas the lowest expression was observed in Dirk with an expression ratio of 0.18. MaxiPak-65 exhibited the highest expression of TaMOR (4.01) while the lowest expression was recorded in Dirk with a 0.06 expression ratio.

Discussion
In this study, we evaluated a set of historical wheat cultivars representing the diversity of Pakistan for changes in RSA over time. Root architectural traits such as lateral root number, volume, root length, density, and surface area aid in water uptake from water-deficient soils [27]. The root system of each of three breeding groups responded differently to water-limited conditions. Pre-1965 cultivars exhibited a decrease in 20 root traits under drought stress as compared to control conditions including length, density, volume, and area related root traits. Whereas cultivars released between 1965 and 2000 showed a positive response to drought in terms of length, volume, and diameter but the number of roots was decreased. The number of roots was increased in post-2000 cultivars under water-limited conditions along with root surface area and diameter. Zhu et al.  [10] analyzed modern high yielding cultivars of northwestern China and showed that primary seminal root length was positively correlated with yield. Roots of higher-yielding modern varieties were simpler,compact and grew deeper than low-yielding modern varieties. Their results were also consistent to this study that the success of wheat breeding for higher yields over past 100 years in northwestern China has been in part due to unconscious group selection on root traits, resulting in smaller, less branched, and deeper roots, suggesting a direction for further increases in crop yield in the future. However, to deliver increases in crop yields, context-dependent optimization of root systems is crucial. For instance, the first cloned gene DRO1 controlling narrow root angle in rice was described to provide yield benefit under limited water conditions [26]. McGrail and McNear [28] also concluded while studying RSA of cultivars representing two centuries of breeding progress that landraces and modern cultivars have contrastingly different RSA, while RSA of intermediate and modern wheat cultivars did not vary significantly.
The response of root growth to water deprivation usually includes growth enhancement of first-and secondorder roots and inhibition of lateral roots growth. When water scarcity is severe, a drought avoidance program is implemented to direct root growth and branching into resource-rich regions [29]. Root biomass in wheat is a multi-trait function including number, length, and diameter of seminal and nodal roots the response of which to drought might be positive, negative, or no response. In our study, the positive response of some important RSA traits of post-1965 cultivars to breeding for drought stress tolerance improved the cultivars ability to transport more assimilates to roots for an efficient root system required for resource uptake. In line with our findings, Ephdaie et al. reported that under severe drought intensity (36%), wheat cv. Pavon76 used a large portion of plant-available water to increase its root biomass thus leaving behind a small amount insufficient for grain filling causing a reduction in grain yield [30]. It has been well established while comparing landraces and modern cultivars that total root biomass significantly differed between Turkish landraces and modern cultivars [31], old and modern Mediterranean wheat [32], old and modern American wheat [33], drought tolerant landraces and CIMMYT-derived wheat [34], and drought tolerant landraces and modern American wheat [31]. All these studies have consensus that landraces or older cultivars have greater biomass compared to modern cultivars.
Root angle is also an important drought-adaptive trait that directs the horizontal and vertical distribution of roots into the soil. A strong link between steep root angle and deep rooting has been reported in wheat [35]. The narrow, compact and deep-rooted architecture of higher-yielding cultivars appears to minimize water use early in the season and subsequently enhance access to water alter in the developmental stages [10]. It has also been reported that introduction of DRO1 into a shallow rooting rice cultivar enabled the cultivar to increase deep-rooting and yield under drought stress conditions [15]. Under water-limited conditions, there was a strong correlation of root depth to other favourable root traits. For instance, root depth was positively correlated to StAF, GD, GL, GY, GPS under drought stress. Zhu et al. [10] also identified negative correlation between seminal root growth angle and grain yield in modern wheat cultivars, which indicated deeper roots likely favored resource acquisition for enhancing grain yield. We also found a positive correlation between depth and StAF in both well-watered and water-limited conditions. Under drought stress, pre-1965 cultivars exhibited a decrease in root depth contrary to other breeding groups after 1965. Progressive enhancement of rooting depth in 1965-2000 and post-2000 cultivars indicates enhanced tolerance to drought stress. Hermanska [36] reported a highly significant correlation between root system size and grain yield. Therefore, the enhanced grain yield in post-2000 cultivars might be attributed to improved root systems tailored to extract water and nutrients from deep soil layers under water-limited conditions. Root growth angle increases towards gravity due to greater expression of DRO1. This gene enabled the scientists to cope with the drought problems by allowing the roots to penetrate deeper thus aiding in yield enhancement even under limited water supply [14]. In our gene expression studies, we found a strong expression of the DRO1 gene in all cultivars tested that supports our findings of increased depth under drought stress as compared to control in cultivars released after 1965.
Although extensive efforts have been made to identify the association of functional genes with the aboveground phenological traits in wheat, and functional markers have been used by various groups to tag favorable alleles in germplasm [22]. Although various Genomewide association studies are available for genome-wide association of SNP markers with RSA. However, less efforts have been made to identify the useful allelic variation with RSA, to further use such information for breeding desirable RSA. RSA may hold the key for the "second green revolution" therefore the effect of semi-dwarfing alleles (Rht) on the root system is of prime importance [37]. Generally, the effects of Rht1 alleles on root systems are less clear, with studies in different growing conditions producing contradictory results. We found a significant association of reduced height alleles (Rht) with some root traits. Rht8 allele is linked to reduced plant height and the highest frequency (90.9%) was recorded in post-2000 cultivars. The results are consistent with our field data where the least plant height was recorded in post-2000 cultivars. The gibberellin (GA-3) insensitive green revolution allele Rht-D1b linked to dwarfism was also frequent (87.8%) in post-2000 cultivars. Similar to our findings, Hurd [38] reported that semi-dwarfing wheat lines had larger root systems as compared to tall control. Whereas contradictory to Laperche et al. [39], we found no significant impact of Rht-B1 on root architecture that might be attributed to differing growing conditions [40].
The low heritability of RSA components traits a major challenge because most of the traits are controlled by many genes with minor effects [41]. Very few QTL with large effect are known for RSA in wheat, of which a QTL on chr5B, qRDM-5B, is prominent and is suggested to be used for selection of breeding germplasm [16]. This QTL was significantly associated with dry root weight (dry root biomass) under both control and drought stress conditions in spring wheat cultivars of Pakistan. The favorable haplotype allele frequency was higher in cultivars released in rainfed areas, however, some old cultivars like C-217 and C-518 also carried favored allele. Wheat LBD (LARGE ORGAN BOUNDARIES DOMAIN) gene TaMOR also plays a significant role in wheat root development and improvement of root systems [42]. Of the 12 genotypes tested, we found the highest RNA-seq based TaMOR expression in five wheat cultivars released in 1965 of which Zincol-2016, GA-2002, and Chakwal-50 belong to the post-2000 breeding group. In line with our results, Li et al. [43] reported that overexpression of TaMOR in rice plants contributed to the larger root system and higher grain yield in rice plants.

Conclusion
It is concluded from our findings that the historical bread wheat cultivars released before the green revolution (1965) have poor root systems as compared to cultivars released thereafter. Cultivars of the post-1965 group have improved root systems under both well-watered and water-limited conditions and thus can tolerate climatic fluctuations. A significant improvement in droughtadaptive traits such as the depth and steep angle frequency was prominent feature of modern cultivars. The favorable allele of a major root biomass QTL, qDRM-5B, was present in low frequency, and its positive selection in breeding could improve drought adaptive RSA. This progressive improvement in RSA is also linked to green revolution Rht genes along with other phenology-related genes such as TaLTPs, TaSus-2B, TaMOR, and DRO1 that have significant allelic effects on RSA traits. Furthermore, correlation analysis of root and yield-related indicates a strong association of root systems with grain yield. Thus, we suggest that it is crucial to integrate the knowledge of promising root systems and root traits into breeding programs to develop climate-resilient varieties.
Additional file 1: Table S1. The raw data for RT-PCR experiment.